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Abstract. We discuss charm production at the LHC. The production of single cc pairs 
is calculated in the ^-factorization approach. We use Kimber-Martin-Ryskin unintegrated 
gluon distributions in the proton. The hadronization is included with the help of Peter- 
^jqJ son fragmentation functions. Transverse momentum and pseudorapidity distributions of 

r—{ _ charmed mesons are presented and compared to recent results of the ALICE, LHCb and 

ATLAS collaborations. Furthermore we discuss production of two pairs of cc within a 
simple formalism of double-parton scattering (DPS). Surprisingly large cross sections, 
comparable to single-parton scattering (SPS), are predicted for LHC energies. We discuss 
perspectives how to identify the double scattering contribution. We predict much larger 
cross section for large rapidity distance between charm quarks from different hard parton 
scatterings compared to single scattering. 

Cl[ 1 Transverse momentum spectra of open charm mesons at LHC 

D . 

Recently ALICE, LHCb and ATLAS collaborations have measured inclusive transverse momentum 
spectra of open charm mesons in proton-proton collisions at yfs — 7 TeV II 1 12 13 1 . These measurements 
are very interesting from the theoretical point of view because of the collision energy never achieved 
before and unique rapidity acceptance of the detectors. Especially, results from forward rapidity region 
2 < y < 4, obtained by the LHCb as well as ATLAS data from wide pseudorapidity range \tj\ < 2.1 
can improve our understanding of pQCD production of heavy quarks. 

The inclusive production of heavy quark/antiquark pairs can be calculated in the framework of the 
fc,-factorization [4 1 . In this approach transverse momenta of initial partons are included and emission of 
gluons is encoded in the so-called unintegrated gluon, in general parton, distributions (UGDFs). In the 
leading-order approximation (LO) within the £,-factorization approach the differential cross section 
for the QQ can be written as: 

dcr v r d 2 k\ , d 2 ki , 1 



dy l dy 2 d 2 p l ,d 2 p2t 



yj J n n \6n l (xiX2sY 
x 5 2 (k u + k 2 , f - pi, f - p 2 , f ) Ti(xi,k 2 t ) Tj(x 2 , k 2 2l ) , 

where Ti(x\ , k 2 ) and < Fj(x2, k\ t ) are the unintegrated gluon (parton) distribution functions. 

There are two types of the LO 2 — > 2 subprocesses which contribute to heavy quarks production, 
99 — * QQ an d qq ~^ QQ- The first mechanism dominates at large energies and the second one 
near the threshold. Only gg — > QQ mechanism is included here. We use off-shell matrix elements 
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corresponding to off-shell kinematics so hard amplitude depends on transverse momenta (virtualities 
of initial gluons). In the case of charm production at very high energies, especially at forward rapidities, 
rather small jc-values become relevant. Taken wide range of x necessary for the calculation we follow 
the Kimber-Martin-Ryskin (KMR) [5 | prescription for unintegrated gluon distributions. More details 
of theoretical model can be found in Ref. J6). 

The hadronization of heavy quarks is usually done with the help of fragmentation functions. The 
inclusive distributions of hadrons can be obtained through a convolution of inclusive distributions of 
heavy quarks/antiquarks and Q — ♦ h fragmentation functions: 

d(r{y h ,p t , h ) _ r l dz D d <^ Q (yQ'P<.<2) 
dy h d 2 p,j, Jo z 2 dy Q d 2 p uQ 

Pl.Q=Pl.lil~ 

where p t Q = — , where z is the fraction of longitudinal momentum of heavy quark carried by meson. 
We have made approximation assuming that yQ is unchanged in the fragmentation process. 

In Fig.Q]we present our predictions for differential distributions in transverse momentum of open 
charm mesons together with the ALICE (left panel) and LHCb (right panel) experimental data. The 
uncertainties are obtained by changing charm quark mass m c = 1.5 ± 0.3 GeV and by varying renor- 
malization and factorization scales p 2 = gm 2 , where f e (0.5; 2). The gray shaded bands represent 
these both sources of uncertainties summed in quadrature. Using KMR model of UGDFs we get very 
good description of the experimental data, in both ALICE and LHCb cases. Here, we also compare 
central values of our LO A:, -factorization calculations (solid line) with NLO parton model (dashed 
line) and FONLL |7 1 predicitons (long-dashed line). All of these three models are consistent and give 
very similar results. The only difference is obtained at very small meson p,'s (below 2 GeV), where 
transverse momenta of initial gluons play a very improtant role. 

In Fig |2] we show transverse momentum (left panel) and pseudorapidity (right panel) spectra of D ± 
mesons measured by ATLAS. The representation of theoretical results and uncertainties is the same 
as in Fig. Q] In contrast to the ALICE midrapidity measurements, here the experimental data points 
can be described only by the upper limit of our theoretical predictions. Therefore one can conclude, 
that covering wider range of (pseudo)rapidities (getting larger rapidity differences between produced 
quark and antiquark) the theoretical description of measured data becomes somewhat worse. 



(2) 




Fig. 1. Transverse momentum distributions of D ± and mesons together with the ALICE (left) and LHCb 
(right) data at ->/s = 7 TeV. Predictions of LO A:, -factorization together with uncertainties are compared with NLO 
parton model and FONLL calculations. 
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Fig. 2. Transverse momentum (left) and pseudorapidity (right) distributions of D ± mesons with the ATLAS 
experimental data at yfs = 7 TeV. Predictions of LO ^-factorization together with uncertainties are compared 
with NLO paton model and FONLL calculations. 



2 Double charm production via Double Parton Scattering 



The mechanism of double-parton scattering (DPS) production of two pairs of heavy quark and heavy 
antiquark is shown in Fig. [3] together with corresponding mechanism of single-scattering production. 
The double-parton scattering has been recognized and discussed already in seventies and eighties. The 
activity stopped when it was realized that their contribution at center-of-mass energies available then 
was negligible. Nowadays, the theory of the double-parton scattering is quickly developing (see e.g. 
B8I9I10I1 ) which is partly driven by new results from the LHC. 

The double-parton scattering formalism in the simplest form assumes two single-parton scatter- 
ings. Then in a simple probabilistic picture the cross section for double-parton scattering can be written 
as: 

a° PS {pp -> ccccX) = —^—(r sps (pp ccXi) ■ o- SPS (pp ccX 2 ). (3) 

2cr eff 

This formula assumes that the two subprocesses are not correlated and do not interfere. At low energies 
one has to include parton momentum conservation i.e. extra limitations: x\ + xj, < 1 and x 2 + x\ < 
1, where x\ and X3 are longitudinal momentum fractions of gluons emitted from one proton and x 2 
and X4 their counterparts for gluons emitted from the second proton. The "second" emission must 
take into account that some momentum was used up in the "first" parton collision. This effect is 
important at large quark or antiquark rapidities. Experimental data IfTTI provide an estimate of <r e ff in 
the denominator of formula (0. In our analysis we take cr e ff = 15 mb. 




Fig. 3. SPS (left) and DPS (right) mechanisms of (cc)(cc) production. 
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A more general formula for the cross section can be written formally in terms of double-parton 
distributions (dPDFs), e.g. F gg , F qq , etc. In the case of heavy quark production at high energies: 

dCT DPS = — Fgg(XU Jt 3 ,^J, n\)Fgg{X 2 , X^,fl\,n\) X 

lcr eff 

X2,n\)do-gg^ c - c {xi,X4,n\) dx\dx2<ix^dx A . (4) 

It is physically motivated to write the dPDFs rather in the impact parameter space F gg (x\, x%, b) = 
g(xi)g{x2)F(b), where g are usual conventional parton distributions and F(b) is an overlap of the matter 
distribution in the transverse plane where b is a distance between both gluons [ 12 1. The effective cross 
section in (01 is then l/o" e // = J d 2 bF 2 (b) and in this approximation is energy independent. 

In the left panel of Fig.|4]we compare cross sections for the single cc pair production as well as 
for single-parton and double-parton scattering cccc production as a function of proton-proton center- 
of-mass energy. At low energies the single cc pair production cross section is much larger. The cross 
section for SPS production of cccc system [13] is more than two orders of magnitude smaller than 
that for single cc production. For reference we show the proton-proton total cross section as a function 
of energy as parametrized in Ref. [14|. At low energy the cc or cccc cross sections are much smaller 
than the total cross section. At higher energies the contributions approach the total cross section. This 
shows that inclusion of unitarity effect and/or saturation of parton distributions may be necessary. At 
LHC energies the cross section for both terms becomes comparable. This is a new situation when the 
DPS gives a huge contribution to inclusive charm production. 




Fig. 4. Total LO cross section for single cc pair and SPS and DPS cccc production as a function of center-of-mass 
energy (left panel) and differential distribution in rapidity difference (right panel) between c and c quarks at ^[s = 
7 TeV. Cross section for DPS should be multiplied in addition by a factor 2 in the case when all c (c) are counted. 
We show in addition a parametrization of the total cross section in the left panel. 

In the right panel of Fig. [4] we present distribution in the difference of c and c rapidities ydtff = 
He - He- We show both terms: when cc are emitted in the same parton scattering (ciq or C3C4) and when 
they are emitted from different parton scatterings {c\cn or C2C3). In the latter case we observe a long 
tail for large rapidity difference as well as at large invariant masses of cc. 

In particular, cc (or cc) should be predominantly produced from two different parton scatterings 
which opens a possibility to study the double scattering processes. A good signature of the cccc final 
state is a production of two mesons, both containing c quark or two mesons both containing c antiquark 
(D°D° or/and D°D°) in one physical event. A more detailed discussion of the DPS charm production 
can be found in our original paper Ref. fT5l . 

In the present approach we have calculated cross section in a simple collinear leading-order ap- 
proach. A better approximation would be to include multiple gluon emissions. This can be done e.g. 
in soft gluon resummation or in ^-factorization approach. This will be discussed in detail elsewhere 

ma. 
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